Generation of a spin-squeezed state with trapped ions using a dressing field 
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We propose a method for generating a spin-squeezed state that is a symmetric Dicke state, with 
trapped ions using only global access. The eigenstates of the ions under a strong dressing field 
become symmetric Dicke states and the M0lmer-S0rensen interaction selectively couples one of 
them to an initially populated auxiliary state. A | DJn) state, which is maximally spin squeezed, 
can be generated with high fidelity using only square pulses. Using an adiabatic technique, the ideal 
maximally spin-squeezed state is generated. 

PACS numbers: 



The generation of genuine multipartite entanglement 
is of central importance for quantum information sci- 
ence and quantum information processing. Entangled 
states are resources for quantum computation and can 
also be used to improve the sensitivity beyond the 
shot noise limit of an atomic clock or in other precise 
measurements [lj. In particular, Greenberger-Horne- 
Zeilinger states and the symmetric Dicke states have 
been extensively investigated [J 0]. The symmetric Dicke 
states are given by 



V nCn 



where Pk are permutation operators and „C m = 
n\/[m\(N — to)!]. Among the symmetric Dicke states, 
the maximally spin-squeezed state (l-D^n)) has the high- 
est sensitivity forprecise measurement and achieves the 
Heisenberg limit [2|. High-sensitivity measurements us- 
ing entangled states have been demonstrated by a single 
pair [1] and two_pairs[3, HJ of photons and a single pair 
of trapped ions[6[. Furthermore, precise measurement 
using the large-scale entanglement of twin matter waves 
has recently been reported^, H|, though its sensitivity 
has not reached the Heisenberg limit. The generation 
of large-scale entanglement with high fidelity remains a 
significant challenge. 

Currently, large-scale entanglement is generated exper- 
imentally by using a common phonon mode with trapped 
ions. In this method, a phonon mode is excited with ei- 
ther an array of laser pulses irradiated locally [ij [l(| o r by 
off- resonant global irradiation (geometric phase gate 11 
and M0lmer-S0rensen gate[l2j|)- Using local pulse arrays, 
a W state (| of eight ions is generated [10|. Using off- 
resonant lasers, the generation of a Greenberger-Horne— 
Zeilinger state of 14 ions has been demonstrated (l3j| . As 
the number of atoms used is increased, local access be- 
comes more difficult and complicated. It is then impor- 
tant that the generation of entangled states does not re- 
quire local access. 

The generation of a symmetric Dicke state was 
proposed [H| and demonstrated [HJ using an adiabatic 
technique. However, this proposal requires a phononic 



Fock state for the initial state. Also, to generate the max- 
imally spin-squeezed state, a large Fock state is needed, 
which is a challenge for achieving high fidelity. 



In this paper, we propose a method for generating a 
maximally spin-squeezed state using only global irradia- 
tion by a laser, radio-frequency or microwave field. The 
combination of a dressing field and the M0lmer-S0rensen 
interaction [3] has been reported experimentally fl7|| for 
two ions. This work presents a theoretical extension to 
a large number of ions. The method can generate the 
large-scale symmetric Dicke state, especially the maxi- 
mally spin-squeezed state (l-D^n))- 



We first consider the situation depicted in fig. 1(a). 
2N ions are trapped and cooled to the ground state in 
a linear RF trap whose normal mode frequency of the 
center of mass mode is uj rn . There are three long-lived 
states, namely an ancillary state, |a), and two qubit 
states, {||), |t)}- The interaction between these states 
can be expressed by the following Hamiltonian in the in- 
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FIG. 1: (a) Energy-level diagram of trapped 2N ions. oj a is 
the frequency of the transition |o) o |4-)> w m is the normal 
mode frequency and S r ,b is the amount of detuning from the 
sideband transition, (b) Whole Hamiltonian approximated as 
a single multi-level ladder. 
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teraction picture, 
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where a and are the annihilation and creation oper- 



ators of the motional mode, respectively, and <j\ and 



a[ l) and a 2 l > ) are the spin flip operators between 

the {||), |t)} ({|a), ||)}) states of the i-th ion. The first 
term (Hi) gives rise to dressed states of the qubit, ex- 
pressed as |+> = (|f) + ||))/V2 and |-) = (|t) - \i))/V2. 
Oi is the Rabi frequency of the transition between the 
qubit states. The second term in the Hamiltonian (H2) 
corresponds to an MS type interaction 17] ■ 77^2 is the 
Rabi frequency of the sideband transitions and 5 r j, (S r = 
w„- (%,- u>i r , 5 b = con, - u a - io m where uj tr ,ib are the 
frequency of lasers) is the amount of detuning from the 
sideband transitions. To make the MS interaction an 
effective Ising interaction, the detunings 5 r .b are much 
larger than the Rabi frequencies 77^2 and the effective 
Ising interaction is at resonance when 5 r = Sb = S\Vl\. 
We introduce the global spin operators S a — J2i &a i 

and J a = J2t ^i,2) where <7^ (a = {x, y, z}) is the pro- 
jection along the a axis of the Pauli matrices. Using this 
notation, the Hamiltonian is rewritten such that 
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The global spin operator S x commutes with the operator 

a 2 

S (= S x + Sy + Si) and its eigen-states, \S, m), namely 
the dressed states, can be expressed by the eigenvalues of 

S and S x . Because J x and S x are global spin operators 
and the Hamiltonian commutes with the permutation op- 
erators of ions Pk, the symmetry of the whole system 
with regard to permutation of ions is preserved during 
the time evolution with the Hamiltonian. We consider 
an initial state \ aa . . . a). Then, the unique completely 
symmetric state, which does not change with arbitrary 
exchanges of ions, can couple to this initial state due to 
the symmetry of the initial state. We define the state 
such that 

\N a ,m) = P sym [\a,a, . . .,a)(g> \ S max ,m)], 

where iV a is the number of ions in the | a) state and S max 
(—(2N—N a )/2) is the maximum value of S when N a ions 
are in the | a) state and P sym is the symmetrizing opera- 
tor, which is the normalized summation of the operators 
proportional to the permutations. Note that the \N a , m) 
states include the symmetric Dicke states (| D 2 g a * +m )) 
with 2N - N n ions. 



We expand the Hamiltonian with these \N a , m) states: 



H, = 



Y, \N a ,m){N a ,m\ 



N a ,m 
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The first term is diagonalized with these states and we ex- 
amine the coupling of the states due to the second term. 
The initial state (|27V, 0)) is not affected by the dress- 
ing field and hence this initial state is an eigenstate of 
the first term of the Hamiltonian with eigenvalue 0. If 
ill » < -''"|- > is satisfied, the initial state can be coupled 
only to the | N a , 0) states by the Hamiltonian because the 
other states have large energy due to the dressing field 
and are greatly detuned from the resonance of the effec- 
tive Ising interaction. The error in this approximation 
is of the order of [^^—/fli] 2 . Then, the whole Hamil- 
tonian can be factorized into a subspace D spanned by 
D = Span{| N a , 0) | N a = 0, 1, . . . , 27V}, which includes 
the initial state and the maximally spin-squeezed state 
1 0, 0) , and the orthogonal complement D^~: 
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We focus on the subspace D and calculate the matrix 
elements of this effective Hamiltonian. Because J 2 is 
constructed by the products of two spin flip operators, 
the matrix elements of Ho are 



= H 



{N a ,0\H D \N' a ,Q) 

Vo (N a = N' a ,N' a ±2) 
= (N a ^N' a ,N' a ±2). 



Then, the interaction on D can be expressed as a single 
ladder (fig. 1(b)). The dimension of the effective Hamil- 
tonian reduces to N+l. 

To calculate exact values of the matrix elements, we 
expand the maximally spin-squeezed state ISmaxjO) with 
the symmetric Dicke states along the z axis (| S, S z = 
m z )): 



I'S'maxj 0) — 



P\S maX) m z J I S'maxi S z — m z)j 



Using this expansion, the Hamiltonian is expressed such 
that 



HN a ,N a 

Smax 
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FIG. 2: (a) Eigenvalues of the Hamiltonian as functions 
of the detuning of the effective Ising interaction with 16 ions. 
In the inset, the solid and dashed curves shows the energy 
difference between the two highest and lowest energy states, 
respectively, which is related to the adiabatic condition, (b) 
Generation of the maximally spin-squeezed state by an adia- 
batic process using 16 ions. The calculation parameters are 
as follows: sweep time: 2 ms, detuning linear sweep: from 
— 2iv x 28kHz to 2tv x +28 kHz, peak strength of the effective 



Ising interaction: 



2S 



2-7T x 3 kHz, and amplitude sweep: 



Gaussian envelope with e 1 full width of 1.3 ms. 
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where V^u are the matrix elements of the effective Ising 
interaction, 

Vff = (J = n/2,j z = l\J x J x \J = n/2,j z = k). 

We numerically analyze the time evolution using this 
Hamiltonian. To generate the large-scale maximally spin- 
squeezed states more effectively, an adiabatic process is 
useful, as in similar multi-level systems (fig. 1(b)) [l^, 18- 
[2lj . This corresponds to the rapid adiabatic passage 
(RAP) technique[22] using the effective Ising interac- 
tion. RAP is a robust population-transfer method using 
a pulse with a time dependent envelope and frequency. 
Fig. 2(a) shows numerically calculated eigenvalues of the 



Hamiltonian Hp against the detuning of the effective 
Ising interaction (A =| S r \ — | 5b |)). We prepared 
the | 2N, 0) state, which corresponds to the curve at the 
top of fig. 2(a), and swept the detuning. The scalable 
maximally spin squeezed state | 0,0) (=| D^)) can be 
ideally generated through RAP. The numerical result is 
shown in fig. 2(b) with 16 ions and the parameters used 
in the calculation are given in the caption. 

In ordinary two-level RAP we can use either of two di- 
rections for the detuning sweep and both directions are 
equivalent with regard to population transfer. However, 
we here use eigenstates of a non-linear Hamiltonian and 
hence the two sweep directions are not equivalent. At 
the resonance condition (A = 0), the energy difference 
between the two lowest energy states and that between 
the two highest energy states are not equal, so that the 
situation is different depending on the direction of the de- 
tuning sweep and on the state that is used (the highest 
energy state or the lowest energy state). The inset of fig. 
2(a) shows the energy difference between the two lowest 
energy states and two highest energy states of fig. 2(a). 
The conditions of adiabaticity differ depending on the di- 
rection of the detuning sweep due to the non-linearity of 
the Hamiltonian and more efficient population transfer 
is achieved with the highest energy state. Generally, the 
adiabatic condition is expressed as A a d t m where 
A a d is the minimal energy separation to the other states 
and t is the typical time for the detuning and amplitude 
sweep. For the 16-ion case (fig. 2) the minimal energy 
separation from the highest energy state is calculated nu- 
merically to be A ac i = 8.02(77f2 2m ax) 2 /25. 

To date, RAP has not been demonstrated for the effec- 
tive Ising interaction, and there are a few technical chal- 
lenges in applying RAP to the effective Ising interaction. 
First, because the effective Ising interaction arises from 
a two photon process [2l|, the strength of the interaction 
is limited to a few kHz in a conventional RF Paul trap 
and a long interaction time is required for RAP. Such 
a long interaction time induces a large decoherence due 
to magnetic fluctuations and other disturbances. Second, 
to form the maximally spin-squeezed states there are two 
limiting conditions regarding detuning: fii 3> ^^f 1 and 
5 r ,b 3> r 1^2- During the detuning sweep these conditions 
must be satisfied at all time. 

Although we can generate multipartite maximally 
spin-squeezed states through an adiabatic process such 
as RAP, it is technically difficult to use the effective Ising 
interaction in an adiabatic process. Here we analyze the 
generation with a single detuning and square pulse. For 
the two-ion case, the Hamiltonian is expressed in the ba- 
sisof {| 0,0), 1 2,0)} as 

h(j]n 2 ) 2 ( 2 V2\ 
,V2 



Hr 
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Note that the diagonal elements are different to each 
other and a shift appears which non-linearly depends on 
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FIG. 3: Time evolution with square pulses, (a)-(c) Evolution 
of the state populations on resonance: (a) 2 ions, (b) 6 ions, 
(c) 20 ions, (d)-(f ) Evolution with optimized detuning of the 
effective Ising interaction: (d) 2 ions, (e) 6 ions, (f) 20 ions. 
Each solid and dashed curve expresses the population of the 
| Mon,0) states and the | 0,0) states, respectively. 
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FIG. 4: Fidelity vs. the number of ions on a log-log plot. A 
high fidelity (> 0.8) can be achieved even for 300 ions. The 
decay of the fidelity is approximately calculated as Fidelity ~ 
0.1013-/V;~°' 04 by fitting these numerical results. 



the number of ions in the | a) states. For the two-ion 
case, we can cancel this shift by adjusting the detuning 
of the effective Ising interaction A, and the ideal maxi- 
mally spin-squeeze state can be generated. The time evo- 
lution with this Hamiltonian is depicted in fig. 3(a) and 
(d) without and with optimized detuning, respectively. 
In other cases, although the non-linear energy shifts are 
never canceled and become complicated, the population 
of maximally spin-squeezed state can be improved with 
optimized detuning [fig. 3(b), (c), (e), (f)] of the square 
pulse. 

We numerically determined the optimized detuning 



[A op t(-/V; n)] that maximizes the fidelity, depending on 
the ion number: 

A opt (iV ion )- 0.225^^^- 7 . 

The maximal fidelity with varying number of ions with 
optimized detuning of the square pulse is shown in fig. 4. 
By fitting the numerical results, the decrease in fidelity 
is found to be in proportion to the —0.04-th power of the 
number of ions and a fidelity of more than 0.9 is achieved 
for 16 ions. 

In summary, we propose and analyze the generation 
of the maximally spin-squeezed state | D^), which is a 
symmetric Dicke state and can be used to achieve Heisen- 
berg limited precise measurements, using only trapped 
ions with global access of the M0lmer-S0rensen (effective 
Ising) interaction and the dressing field. Under the con- 
dition that the dressing field is much stronger than the 
effective Ising interaction, the whole Hamiltonian can be 
divided into a small subspace and we can calculate the 
time evolution of the system for a large number of ions. 
The fidelity of this state can be made unity by introduc- 
ing an adiabatic process. The ideal large-scale maximally 
spin squeezed state can be generated by sweeping the 
detuning of the effective Ising interaction. The effective 
Hamiltonian includes non-linear energy shifts. However, 
when these shifts are canceled by the detuning of the 
effective Ising interaction, the maximally spin-squeezed 
state can be generated with high fidelity, which decreases 
in proportion to the —0.04-th power of the number of 
ions, even using square pulses. 
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